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Magnetic Resonance Force Microscopy can enable nanoscale imaging of spin lifetime. We report 
temperature dependent measurements of the spin correlation time r m of the statistical fluctuations 
of the spin polarization — the spin noise — of ensembles containing ^100 electron spins by this tech- 
nique. Magneto-mechanical relaxation due to spin-cantilever coupling was controlled and spurious 
mechanisms that can affect the spin correlation time of the microscopic signal were characterized. 
These measurements have ramifications for optimizing spin sensitivity, understanding local spin 
dynamics and for nanoscale imaging. 



Magnetic Resonance Force Microscopy (MRFM) [TJ [2] 
can detect magnetic resonance from very small spin en- 
sembles with single electron spin sensitivity [3] . For small 
spin ensembles, statistical fluctuations of the net spin 
polarization P net (7V t - 7V ; )/(7V t + N±) [5 |5] exceed 
the Boltzmann polarization. Spin noise is a topic of in- 
trinsic interest [6] as it reveals fundamental information 
about the microscopic environment around the measured 
spins. Spin relaxation provides a powerful approach to 
probing electronic, magnetic and structural dynamics in 
materials [7 , and plays an important role in Magnetic 
Resonance Imaging (MRI) where Ti- and T2-weighting 
are used to enhance image contrast [8]. The signal-to- 
noise ratio (SNR) — singularly important for high reso- 
lution MRFM imaging — is centrally influenced by spin 
lifetime because it determines the detection bandwidth. 

Here, we report measurements of r m in nanoscale en- 
sembles containing ~100 electron spins. The number of 
resonant spins and the correlation time r m of their fluc- 
tuations are characterized in MRFM experiments by the 
spectral weight and linewidth, respectively. Ideally, the 
spin-lattice relaxation time in the rotating frame, Ti p , 
determines r m [9j [10] . We present systematic measure- 
ments of the evolution of r m with spin modulation depth, 
microwave power and sample temperature. We argue, 
based on these data, that the relaxation time we measure 
in these experiments are due to intrinsic mechanisms. 

Care must be taken to avoid artificially shortening 
the spin correlation time through mechanisms of tech- 
nical origin such as the higher order cantilever oscillation 
modes [9HTT] , violation of adiabaticity [12] [13] , and low 
frequency fluctuations of the field of the micromagnetic 
probe [14]. We avoided these by using mass- loaded can- 
tilevers [15], large cantilever oscillation magnitudes x p k 
and large transverse oscillating magnetic fields H\. We 
find that the temperature dependence of l/r m (T) are in- 
trinsic to the sample and are well explained by phonon 
mediated relaxation processes. 

Our experiments were performed in vacuum between 
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FIG. 1: Power spectral density of the spin noise. The full 
width at half maximum and the area of the fitted Lorentzian 
are 0.30 Hz and 162 mHz 2 respectively. They correspond to 
Tm — 1-06 s and signal energy of 534 aN 2 , which is equivalent 
to 302 electron spins in 80 nm 3 voxel with 1.3 G/nm field 
gradient. Data taken with tip-sample separation d — 350 nm. 



4.2 and 40 K on an optically polished piece of vitreous 
silica (see Ref. 16 for details). We measure electron spins 
present at a density of ~ 6 x 10 1T cm -3 . These spins 
reside in silicon dangling bonds associated with oxygen 
vacancy defects known as E' centers, which are produced 
by 60 Co gamma irradiation [T7H2T] . The sample is ther- 
mally anchored to a temperature controlled copper block. 
The IBM-style ultrasoft cantilever we used has a spring 
constant k ~ 0.1 mN/m and a mass- loaded tip to sup- 
press tip motion [15] arising from thermal excitation of 
higher order cantilever oscillation modes [9]. The probe 
magnet is a SmCos particle glued to the cantilever and 
ion-milled to a tapered end whose size is about 300 x 600 
nm 2 . It has coercivity and anisotropy fields greater than 
1 T at low temperature, thus avoiding spin relaxation 
induced by fluctuations of the probe magnetic field [14]. 
The cantilever frequency f c is 3062.15 Hz with the tip at- 
tached. The transverse oscillating (2.162 GHz) magnetic 
field Hi is generated by a superconducting microwave 
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resonator [2Z . The experiments were performed with no 
external magnetic field applied. 

We used the iOSCAR protocol [4] to excite mag- 
netic resonance and measured the resulting cantilever fre- 
quency shift Sf c resulting from the modulated magnetic 
interaction between the electron spins and the micromag- 
netic probe on the cantilever. Random and uncorrelated 
spin noise leads to a Lorentzian frequency dependence of 
the power spectral density Sf c of these frequency shifts 
(see Fig. [T]) as in the random telegraph signal model 
[3 [23]: 



s fc = 



1 + 4ttV2 



frr 



(1) 



where f m is the iOSCAR modulation frequency and e/ 
is the average frequency-signal energy from N resonant 
electron spins. The area under the Lorentzian in Fig. [I] 
is 162 mHz 2 ; this gives a force signal energy e of 534 
aN 2 ; the two are related by e = (7rkx pk /2f c ) 2 6f [4l [24]. 
The measured tip field gradient is about 1.3 G/nm, so 
the statistical polarization is due to approximately 302 
electron spins (y/~N = 17 A) in a ^ (80 nm) 3 detected 
volume. The noise floor, 13 aN/^/Hz, is primarily due to 
thermal force noise and corresponds to spin sensitivity 
of ~ 100 electrons in a 1 Hz bandwidth. Hereafter both 
r m and e are taken from a fit to the single side-band 
power spectral density obtained by means of a software 
lock-in amplifier with a bank of low pass filters [3] [5] to 
improve SNR. Most of the data points take about 1 hour 
for averaging. 

The correlation time r m is determined by the relax- 
ation time in the rotating frame T\ p averaged over the 
distribution of effective field frequencies experienced 
during an adiabatic inversion cycle [9] [10]. In the absence 
of excess low frequency spin fluctuations, T\ p approaches 
Ti [25j [26] . If the spin spends most of its time far off- 
resonance, that is, if either the microwave frequency mod- 
ulation or the product of the cantilever oscillation ampli- 
tude and the probe gradient is large (ensuring that the 
extremum of the time- varying effective magnetic field in 
the rotating frame i7 e ff is much larger than Hi), and if 
the adiabatic condition is satisfied, then r m in iOSCAR 
should approach T\. 

We explored the dependence of l/r m on x p k at three 
temperatures (see Fig. 2|. Similar to [26 , we find l/r m 
decreases asymptotically to the temperature dependent 
intrinsic relaxation rate: l/r m = /3(T)x^ k + T 1 _1 (T) 
where —1 < a < —0.7 (dashed lines). As the resonant 
slice sweeps through larger volumes with increasing x p k, 
e increases linearly (lower panel). 

Both violation of adiabaticity and magnetic field fluc- 
tuations due to higher order cantilever modes [9UT3] can 
limit r m . To ensure our results are free of such arti- 
facts, we studied the dependence of l/r m on microwave 
power P^ w : Fig. [3] shows l/r m to be independent of P^ w 
for P^ w > 0.4 mW. At low power, l/r m increases due 
to violation of adiabaticity (black dashed line) or other 
mechanisms (the shoulder near 0.15 mW). The measured 
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FIG. 2: Dependence of relaxation rate l/r m on cantilever os- 
cillation amplitude x p k at three different temperatures. l/r m 
decreases asymptotically toward T 1 _1 (T) as x p k increases. 
The dashed lines are the fits to phenomenpower law behavior 
l/rm = P(T)x* k + T' 1 ^) where -1 < a < -0.7. Lower 
panel: the signal energy increases linearly with x p k as ex- 
pected. In this experiment, d — 200 nm, inw — 2.51 mW and 
f m = 45 Hz for T = 18.8 K, 19.9 Hz otherwise. 



signal energy e increases and saturates as P^ w increases. 
Thus we can access a measurement parameter regime in 
which r m measures intrinsic relaxation. 

To understand the relaxation mechanism, we measured 
the temperature dependence of l/r m . These measure- 
ments are presented in Fig. [4] To avoid spurious reduc- 
tion of r m and thus ensure that r m (T) represents Ti(T), 
Xpk and P^ w were kept at 85 nm and 2.51 mW respec- 
tively. The cantilever was thermally isolated from the 
sample and there was no observable increase in the ther- 
mal force noise with change in sample temperature. The 
data are well described by the function 



AT + C 



pTlm/T 



( e T lm /T _ !)2 



(2) 



where A, C and T\ m are fitting parameters. For T < 
16 K, direct phonon absorption or emission dominates 
l/r m , and the linear temperature dependence reflects the 
phonon mode occupancy n = (exp (huj^m/kBT) — 1) _1 ex 
ksT /hj0^ m for T <C hu^ m /kB, where uo^ m is the mi- 
crowave angular frequency, H and ks are the Planck 
and Boltzmann constants respectively. For T > 16 
K, the data for l/r m are well fit by the expression 
C e ^im/T '/( e Ti m /T _ ^2 [isj [19], the signature of two- 
phonon Raman process, in which phonons with frequency 
fim — ^sTim/h are created and annihilated. Unlike the 
usual Raman process in which all available phonon modes 
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FIG. 5: Power law dependence of the signal energy on r" 1 : 
e oc (l/r m ) a plus an offset with the sample temperature (red 
circles) and microwave power (black squares) as implicit pa- 
rameters. Fits to data give a — —0.56 ±0.05 and —0.54 ±0.04 
for varying the sample temperature and microwave power (us- 
ing only l/r m > 0.7 Hz) respectively. 



FIG. 3: Relaxation rate l/r m and signal energy e versus mi- 
crowave power Pftw. At low microwave power, extraneous 
mechanisms such as violation of adiabaticity (black dashed 
line [TQl [27] ) will increase l/r m . The red dotted line shows 
the intrinsic l/r m o phonon mediated relaxation process. The 
measured signal energy e saturates at an intrinsic value eo as 
microwave power increases (green dashed line). 
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FIG. 4: Relaxation rate l/r m versus sample temperature. 
The red dotted line is the best fit to the direct (AT) and Ra- 
man local mode Ce Tlm ^ T / (e Tlm ^ T — l) 2 relaxation processes. 
The dashed lines show the two processes independently. Data 
taken with tip-sample separation d = 200 nm. 



below the Debye frequency can induce electron spin re- 
laxation, this mechanism involves thermal excitations of 
the local mode of the oxygen vacancy defect only. Sil- 
ica with defects induced by neutron-irradiation [28j [29] 
and hydrogen impurities [30[[3T] behave similarly. Acous- 
tic attenuation [32j [33] , infrared and Raman studies also 
support the local mode model. 

Since the local mode frequencies are expected to de- 
pend on the specific type of quartz, the agreement be- 
tween fitted T\ m = 124 ± 18 K and reported values 
[T7HT9] is satisfactory. However, the fitted value A = 
0.112 ± 0.004 Hz is roughly 20 times larger than that 
reported in Ref. [18] The direct process depends on the 
Zeeman splitting, so the discrepancy may arise from the 
different microwave frequencies used. We also find the 
scaling of A deviates from the expected oo 2 w behavior; 
we note similar behavior has been reported [29, 34 and 
suggested to result from a cross relaxation process. 

Fig. [5] shows a remarkable anticorrelation between 
l/r m and e. When e is caused to vary through its de- 
pendence on either P^ w or temperature, we find l/r m 
varies linearly with e a with a ~ — 0.5. A similar depen- 
dence can be found in Ref. [13] We expect e to decrease 
for l/r m ^> f m due to decreased sensitivity to variations 
of the spin magnetization occurring within a single mod- 
ulation period [35]. This cannot explain the observed 
variation of l/r m well below f m (45.1 and 21.5 Hz for 
the P^ w and temperature scans, respectively); the power 
a ~ —0.5 is also not consistent with this origin. Same 
set of filters were used in the cantilever measurement and 
control circuits in these data. 

We have demonstrated the ability to measure local spin 
relaxation times using ultrasensitive MRFM. This points 
to the capability for microscopic measurement of the spa- 
tial variation of spin dynamics. This can provide new 
insight into materials such as superconducting cuprates 
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where intrinsic inhomogeneity plays a central role [36] , 
and could provide essential understanding in technologi- 
cally important phenomena such as spin coherence, spin 
transport and quantum information processing. Further- 
more intrinsic correlation times can provide a mechanism 
for enhancing information content of images through re- 
laxation rate contrast in analogy to T\- and T2-weighted 
Magnetic Resonance Imaging [8]. We find, importantly, 
that care must be taken to understand and account for 
the influence of size sample size on the measured spin 
dynamics. 

Understanding and manipulating r m also has implica- 
tions for MRFM sensitivity. Spin noise detection SNR 
depends on r m [5] because of the trade off between the 
averaging counts and lock-in detection bandwidth. Ref.|5] 
uses 7r/2 rf pulses to randomize the spins and hence re- 
duce r m to the optimal point for maximum SNR. As a 
consequence of the strong field gradient this required very 
broadband rf field which was achieved through trains of rf 
pulses. Our result suggest the optimal r m can be achieved 
by controlling the sample temperature. 



We have studied the dependence of spin relaxation in 
few-electron-spin ensembles on x p k and microwave power 
to understand the effect of time spent near resonance H\ 
amplitude respectively. We have measured the intrinsic 
correlation time r m of the spin noise using the statistical 
polarization signal in ensembles of ~ 100 electron spins in 
vitreous silica. Relaxation is due to coupling of spins to 
phonons through either a direct (single phonon) process 
at low temperature or a Raman process at higher temper- 
atures. This demonstrates the capability for microscopic 
measurement of electron spin dynamics, an important 
quantity for understanding fundamental characteristics 
of electronic systems. Furthermore, understanding and 
controlling r m will be important for future MRFM imag- 
ing applications and sensitivity optimization. 
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